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We study complex time series (spike trains) of online user communication while spreading messages 
about the discovery of the Higgs boson in Twitter. We focus on online social interactions among 
users such as retweet, mention, and reply, and construct different types of active (performing an 
action) and passive (receiving an action) spike trains for each user. The spike trains are analyzed by 
means of local variation, to quantify the temporal behavior of active and passive users, as a function 
of their activity and popularity. We show that the active spike trains are bursty, independently of 
their activation frequency. For passive spike trains, in contrast, the local variation of popular users 
presents uncorrelated (Poisson random) dynamics. We further characterize the correlations of the 
local variation in different interactions. We obtain high values of correlation, and thus consistent 
temporal behavior, between retweets and mentions, but only for popular users, indicating that 
creating online attention suggests an alignment in the dynamics of the two interactions. 
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I. INTRODUCTION 


In recent years, online social media (OSM) have be¬ 
come a major communication channel, allowing users 
to share information in their social and professional cir¬ 
cles, to discover relevant information pre-filtered by other 
users, and to chat with their acquaintances. In addition 
to their practical use for individuals, OSM have the ad¬ 
vantage of generating a rich data set on collective social 
dynamics, as social relations among individuals, tempo¬ 
ral properties of their interactions, and their contents are 
automatically stored. The study of these digital foot¬ 
prints has led to the emergence of computational social 
science, allowing to quantify at large-scales our political 
ideas and preferences [I|, to discover roles in social net¬ 
work ii, to predict our health Q and personality [3, 
and to determine external effects on online behavior g. 
Importantly, in OSM, users are at the same time both 
actors and receivers and therefore the amplification of 
a trend originates from the interplay between influenc¬ 
ing [3,0 and being influenced [ol-fl^ . 

A crucial aspect of OSM and more generally of hu¬ 
man behavior is the underlying complex dynamics [H- 
[i3- The time series of user activities, e.g. posting a 
tweet and replying to a message, are quite distinct from 
uncorrelated (Poisson random) dynamics in the pres ence 
of burstiness temporal correlations 12iL 2 ^. 

and non-stationarity of human daily rhythm [231 |24| . 
which has significant implications. Diffusion on a tem¬ 
poral network cannot be accurately described by models 
on static networks and consequently the process presents 
non-Markovian features with strong influence on the time 
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required to explore the system [^[^. Furthermore, the 
dynamics drives a strong heterogeneity observed in user 
activity [13, and user/content popularity [29l - [3ll| . 
Specifically, in Twitter, the heterogeneity in popularity 
has been observed and quantified in different ways by the 
size of retweet cascades, i.e. users re-transfer messages to 
their own followers with or without modifying them (33 - 
ll^l or by the number of mentions of a user name, identi- 
fied by the symbol @, in other people’s tweets [37| . 

In this paper, we focus on the dynamics of social in¬ 
teractions taking place when diffusing rumors about the 
discovery of the Higgs boson on July 2012 in Twitter [s^. 
Our main goal is to find connections between the sta¬ 
tistical properties of user time series established on the 
same subject, e.g. the announcement of the discovery 
of the Higgs boson, and their activity and popularity. 
To this end, we analyze tweets including social interac¬ 
tions, such as retweets of a message (RT), mentions of a 
user name (@), and replies to a message (RE). For each 
type of the interactions, a user can either play an ac¬ 
tive, e.g. retweeting, or a passive, e.g. being retweeted, 
role. Therefore, we characterize each user by 8 time se¬ 
ries: one active and one passive time series for each of 
the 3 types of interaction as well as for the aggregation 
of all interactions, as illustrated in Fig. |T] Active time 
series are denoted as WHO and passive time series are 
defined by WHOM. We then investigate whether the sta¬ 
tistical properties of each signal is a good predictor for 
the activity and popularity of a user. 

The following sections are organized as follows. In sec¬ 
tion im we describe the data set and provide basic sta¬ 
tistical properties of who and whom time series. In sec- 
tion lHIl we introduce a technique dedicated to the analy¬ 
sis of non-stationary time series, so-called local variation, 
originally established for neuron spike trains [33 - l43| and 
recently has been applied to hashtag spike trains in Twit- 
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ter [ 4 ^, . In section lIVi we search for statistical rela¬ 

tions between local variation and measures of popularity 
of a user. Finally, section fVl summarizes the key results 
and raises open questions. 


WHO 



WHOM 



Pjj\ popularity of users 


Figure 1. Illustration of communication in Twitter. Users in 
who interact in time with users in whom by retweeting (RT) 
the messages and mentioning (@) the user names of whom 
in a message and replying (RE) to the messages from whom. 
Quantifying temporal patterns in time series of who users with 
various ranges of the activity of users au and of whom users 
by increasing the popularity of users pu is the main scope of 
this paper. 


II. ACTIVITY AND POPULARITY OF USERS 


Data Set. As a test bed, we consider the publicly 
available Higgs Twitter data set [H, , first collected 

to track the spread of the rumor on the discovery of the 
Higgs boson via RT, @ or RE. The data set is composed of 
tweets containing one of the following keywords or hash- 
tags related to the discovery of the Higgs boson, “Ihc”, 
“cern”, “boson”, and “higgs”. The start date is the 1st 
July 2012, 00:00 am and the final date is the 7th July 
2012, 11:59 pm, which covers the announcement date of 
the discovery, the 4th July 2012, 08:00 am. All dates and 
timestamps in the data are converted to the Greenwich 
mean time. Detailed information on the data collection 
procedure and basic statistics can be found in Ref. [s^. 

In total, the data is composed of 456,631 users (nodes) 
and 563,069 interactions. Among those, we detect 
354,930 RT, 171,237 @, and 36,902 RE, which shows 
that RT is more popular than the other communication 
channels. For RT interactions, we find 228,560 who and 
41,400 whom users. These numbers are smaller for @, e.g. 
102,802 who and 31,477 whom users, and even smaller 
for RE, with 27,227 who and 18,578 whom users. In 
each case, whom is much lower than who, as expected 
because a small number of users tend to attract a large 
fraction of attention in both friendship and online 

social networks. This observation is confirmed in 

Fig.l where we present the Zipf plots associated to each 
interaction, clearly showing a strong heterogeneity in the 
system. For who users, the frequency of the user com¬ 
munication fu ranks how active users are and measures 
the activity of users ap, on the other hand, for whom 
users, fu quantifies how often the users or their tweets 
are addressed and so gives the popularity of users pu ■ 


Our aim is to examine the dynamics of user communi¬ 
cation in Twitter. We investigate how frequently Twit¬ 
ter users talk to each other on a certain topic, e.g. the 
discovery of the Higgs boson, and identify how complex 
dynamic patterns of the communication evolve in time. 
To this end, we focus on the 3 different types of interac¬ 
tion between users, retweet (RT), mention (@), and reply 
(RE). Twitter users can adopt a tweet of someone and 
use it again in their own tweet blog by RT or contact to 
other users directly by typing user names in a message 
called @ or simply RE to any tweets, e.g. regular tweets, 
retweets, and tweets/retweets including @s. Typically, 
@s and REs are associated to personal interactions be¬ 
tween users, whereas RTs are responsible for large-scale 
information diffusion in the social network and present 
cascades. Here, we count all types of interaction as a 
part of complex information diffusion in Twitter. 

Interactions in Twitter are performed between at least 
two users (for instance, a user can mention several other 
users in a single tweet). Each action is directed and char¬ 
acterized by its timestamp. The users performing the 
action play active roles (who users) and the users receiv¬ 
ing their attention play a passive role (whom users). We 
construct active and passive RT, @, and RE spike trains 
for each user. 


III. LOCAL VARIATION OF WHO AND WHOM 

Communication Spike Trains. We extract salient 
temporal patterns of the user communication time series. 
We evaluate each directed interaction (RT, @, and RE) 
of the users in the pool of who with any users in the 
whom class, shown in Fig. [TJ We don’t check whether 
the whom users participate the conversation in a later 
stage and only construct independent time series of the 
who and whom users. The elements of the time series are 
the timestamps of the data providing us the exact 

time in second of the interaction and the user name or ID 
of the corresponding who and whom users. Ordering the 
timestamps from the earliest to the latest, we generate 
spike trains carrying full story of the communication of 
each user. The resultant user communication spike trains 
are grouped in eight: For each who and whom, the spike 
trains of all interactions together (i) and the spike trains 
of filtered timestamps of RT (ii), @ (iii), and RE (iv). 

Local Variation. A standard way of investigating 
the dynamics of human communication is to examine 
the statistics of the inter-event spike intervals such as its 
probability distribution [T^ . short-range memory coeffi¬ 
cient and burstiness parameter flB | or Fano factor. How- 
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Figure 2. How often who users communicate with whom 
users. Zipf plots describe heterogeneities of users in the types 
of Twitter interaction, e.g. (a) all of retweet, RT, mention, 
@, and reply, RE, (b) only RT, (c) only @, and finally (d) 
only RE. The frequency of the communication fu is measured 
in two-fold: The activity of who users (red squares) au and 
the popularity of whom users (blue circles) pu- The a;—axis 
ranks the users ru from high fu to low values. Each plot 
indicates that Twitter users, who users, more likely contact to 
someone observed in smoother decays of who users, however 
only few whom users are addressed and become popular in 
these communications. 


ever, recent works have showed that further detail anal¬ 
ysis is required to resolve temporal correlations [3l|, [s^ , 
bursts flol-l^ . and cascading driven by circadian 
rhythm WM, complex decision-making of individu¬ 
als ana l54j|. and external factors 0 such as the an¬ 
nouncement of discoveries, as considered in the current 
data [s^. 

To uncover the dynamics of the communication spike 
trains elaborately, we apply the local variation Ly origi¬ 
nally defined to characterize non-stationary neuron spike 
trains [sol - ld^ and very recently has been used to analyze 
hashtag spike trains [^, 1^ . Unlike to the memory coef¬ 
ficient and burstiness parameter [I^, Ly provides a local 
temporal measurement, e.g. at of a successive time se¬ 
quence of a spike train ..., r^-i, r^, r^+i, ..., and so 
compares temporal variations with their local rates 


2 ^ V (R-I-1 “ R) + iji - Ti-l) 


where N is the total number of spikes. Eq. [T] also takes 
the form [4l| 


3 (^Ti+l ~ 


2 ^ 


( 2 ) 


Here, Ar^+i = r^+i — quantifying the forward delays 
and /S.Ti = Ti — Ti-i representing the backward waiting 


within one second. Such events are rare and their im pac t 
on the value of Ly has been shown to be limited |4^ . 
Fig. [3] describes the distribution of Ly, P{Ly) of full 
spike trains all together with RT, and RE for the who 
(a, b) and whom (c, d) users. Grouping Ly based on 
the frequency fu, e.g. the activity of the who users au 
and the popularity of the whom users pu, we examine 
the temporal patterns of the trains in different classes of 
au and pu- For the real data in (a, c), in Fig. EJa), Ly 
is always larger than 1 in any values of au, suggesting 
that all who users contact to the whom users in bursty 
communications. However, in Fig. EJc), we observe dis¬ 
tinct behavior of the whom users and bursts present only 
for low Pu- By increasing pu, Ly tv 1 indicating that 
there is no temporal correlations among the who users 
referring the whom users and Ly is slightly smaller than 
1 for the most popular users, indicating a tendency to¬ 
wards regularity in the time series, as also observed for 
the hashtag spike trains . These observations are sig¬ 
nificantly different for artificial spike trains constructed 
by randomly permuting the real full spike train and so 
expected to generate non-stationary Poisson processes. 
Therefore, all distributions are centered around 1 in this 
case, independently of au and pu, as shown in Figs.EKb, 
d). The randomization and obtaining a null set follow 
the same procedure explained in detail in Ref. (d^ . 

Even though Fig. [3] represents P{Ly) of full spike 
trains, i.e. all interactions together, P{Ly) of individ¬ 
ual RT, and RE communication spike trains describes 
very similar temporal behavior for both the who and 
whom users. Fig. 2] summarizes the detail of P{Ly), 
the mean of Ly, p{Ly) with the corresponding standard 
deviations a{Ly) as error bars, comparatively. The re¬ 
sults highlight that to classify the communication tem¬ 
poral patterns neither the position of the users, whether 
active or passive, nor the types of the interaction, but 
the frequency of the communication fu such as au and 
Pu plays a major role. All Figs. HKa-d), we observe three 
regions: Bursts in low fu, logio(/c/) < 2.5, irregular un¬ 
correlated (Poisson random) dynamics in moderate and 
high fu, logio(/£/) ~ 2.5-3, and regular patterns in very 
high fu, logic(/c/) > 3. This conclusion supports the im¬ 
portance of frequency so time parameter overall human 
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Figure 3. Probability density function of the local variation 
Lv, P{Lv) of who (a, b) and whom (c, d) users in various 
ranges of the two communication frequencies, e.g. au and pu¬ 
la., c) describe the results of the real data. When we only ob¬ 
serve bursty communication patters in who users independent 
of the average user activity frequency {au) in (a), significant 
variations in Lv by increasing the average user popularity 
{pu) are clear in (c). The results prove that popular users 
in Twitter are addressed randomly in time and slightly more 
regular patterns observed in the most popular users. On the 
other hand, (b, d) present the statistics of artificially gener¬ 
ated random spikes serving as a null model and all frequency 
ranges give the distributions around 1, as expected for tem¬ 
porarily uncorrelated signals. 


Figure 4. Mean p of the local variation Ly of the user com¬ 
munication spike trains versus the logarithmic average fre¬ 
quency logic (/}. The results of who users are represented 
by red squares and blue circles describe that of whom users. 
Types of the interaction are investigated in detail: (a) All 
communications of retweet, RT, mention, and reply, RE. 
(b) Only RT. (c) Only @. (d) Only RE. Independent of the 
types of the interaction, the frequency of communication, e.g. 
the activity of users au and the popularity of users pu, designs 
overall communication patterns. While low fu gives bursty 
patterns with Ly > 1, moderate fu indicates irregular uncor¬ 
related (Poisson random) signals, e.g. Ly « 1. For all high 
fu, Ly < 1 presenting the regularity of the communications. 
The error bars show the corresponding standard variations. 


behavior 

We now perform more detail comparison in Fig.jSj how 
Lv of different interactions in the same frequency range 
varies from each other. To this end, we calculate the 
standard z—values in two ways. First, to compare Ly of 
the full spike trains with Ly oi only RT and also with 
Ly of only @ spike trains, and Ly, respectively, we 
introduce 


,j, ^ - Lo{Lv) 

z{fu) = - 1 =^ (3) 

Here, k in superscripts labels the interaction, e.g. either 
RT or @. Precisely, Ly is determined based on a filtered 
spike train composed of the user timestamps of either RT 
or @, as already used in Fig. HKb-c). In addition, pf is 
the mean of Ly, also presented in Fig. HKb-c), and po is 
the mean Ly of the full spike train, given in Fig. SKa). 

In Fig.[5l black squares show z—values of RT and black 
circles describe z—values of @. For the who users in 
Fig.EKa) where Ly only presents bursty patterns (orange 
shaded area) and low ajj, we have small z— values proving 
the agreement of the temporal patterns suggested by Ly 
in the same ajj- However, for the whom users in Fig.jSKb) 
where we have rich values of pu compared to the values 
of au, while z—values are small in bursty patterns (low 


Pu, orange area) as also observed in the who users and in 
regular patterns (high pu, yellow area), larger z—@ value 
(the black circle) is calculated in uncorrelated Poisson dy¬ 
namics (moderate pu, purple area). The disagreement of 
Ly with large z—@ indicates that even though Ly w 1 in 
this region the results of @ are quite sensitive in the same 
Pu, which is not observed in z—RT (the black square). 

Furthermore, we repeat the analysis across communi¬ 
cation channels by comparing temporal patterns of RT 
and @ as follows 


z{fu) 




( 4 ) 


The corresponding z—values, z—@RT are presented in 
green diamonds in Fig. [S] Comparing to the previous 
z—RT and z-@, we now obtain lower values for the who 
users [Fig. [Sl[a)] showing a good agreement between RT 
and @ patterns. Moreover, we have very similar trend 
for the whom users [Fig. EKb)] as before in orange and 
yellow areas and large fluctuations are observed only in 
purple area. 
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Figure 5. Detail comparison between the temporal pat¬ 
terns of different communications in each frequency range. 
While a;—axis is the logarithmic average of frequency, e.g. (a) 
logic (ay) for the who users and (b) logic (py) for the whom 
users, y—axis provides the calculation of three 2 —values (i) 
2 —RT, the comparison of Ly of the full spike train with Lv 
of RT, in black squares, (ii) 2 —@, the same with Lv of @, 
presented in black circles, and (iii) 2 —@RT, the comparison 
between Lv of @ and RT, shown in green diamonds. All 
2 —values are consistent with each other such that except mod¬ 
erate frequency range in (b), e.g. 2 —@ and 2 —@RT, we ob¬ 
serve small 2 concluding that the temporal patterns in the 
similar frequency ranges are in a good agreement. Three dis¬ 
tinct regions are colored due to the discovered patterns in 
calculating Lv in Fig. [T] Orange shaded area describes the 
ranges of the bursty patterns (au and low pu), purple area 
is for the irregular uncorrelated -Poisson random- patterns 
(moderate pu), and yellow area covers the regular patterns 

(highpy). 


IV. CORRELATION OF Lv IN USER 
COMMUNICATION HABITS 


In this final section, our interest turns into building 
new measures to quantify how the local variation Ly fluc¬ 
tuates inside different classes of the frequency, /y. What 
extend temporal communication habits of two users in 
the same /y ranges are dependent on each other is the 
first question we address. Second, we examine whether 
the temporal patterns of the interactions are consistent 


with each other for the same users and how the metric 
varies with increasing /y. 

We consider (fu), the Pearson correlation coeffi¬ 
cient of Ly of two different users selected independently 
from the same fu classes 


Nu , , 

S [^Vi ~ ~ ^{Lyj)] 

- (5) 


Nu 


where <j{Ly.) = \ ~ Li^y^]'^- Here, Ly. and Ly. 


y i=i 

are the local variations of user i and j, respectively, /i’s 
are the corresponding mean values, and Nu is the total 
number of users. Moreover, k and fc' represent all per¬ 
mutations among the full, RT, and @ spike trains. Fur¬ 
thermore, {fu) is evaluated for the who and whom 
users, separately. Therefore, i and j are different users, 
but from the same (who/whom) pool and in the same fre¬ 
quency classes of ay and py, as grouped in Fig. [3] Note 
that before performing Eq. O the corresponding Ly’s 
in the same /y class are ordered from the highest to the 
smallest (or vice versa) not to deform {fu) artificially 
due to the random selection. 

Fig. [3] presents the results of {fu) for the who 
users in (a, b) and the whom users in (c, d). Similar 
to 2 —values performed in the previous Section, we sug¬ 
gest three correlation coefficients: Red (left) triangles de¬ 
scribe blue (right) triangles are for and 


black and green diamonds show the values of ’ . The 
average frequency of the users (/y) in the same class is 
similar but not equal and that is why Figs. [HKb, d) are 
plotted with respect to both the mean frequencies of RT 
and @, e.g. the average activity (ay) and popularity (py) 
of RT and @. All correlations are above 0.85 proving the 
high dependency of the communication patterns of the 
users in the same (/y), independent of the types of the 
interaction. 

We now consider Eq. O with imposing the same user 
and repeat the procedure above for the correlation coef¬ 
ficient 


j.kk' 

' i 


{fu) 




k'{^Vi)][^Vi t^{^Vi)] 


( 6 ) 


Fig.[7]summarizes the results of Eq.lHl While Eigs.[7Ka, 
c) are in parallel with that of Eig. |6] with slightly lower 
correlations for @ (blue right triangles), distinct behavior 
is observed in Eigs.[7Kb, d). Low correlations in Eig. [Tib) 
indicate that the same who users present different tempo¬ 
ral behavior in RT and @. On the other hand, Eig. [T^d) 
shows an interesting temporal habit of the whom users. 
Having no remarkable dependency captured in low pop¬ 
ular users, we show that the correlation increases with 
(py) describing that the popular users are addressed in 
RT and @ in a temporarily similar procedure. 
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Figure 6. Linear correlations of Lv of user pairs: The stan¬ 
dard Pearson correlation coefficient quantifies the dependency 
on the temporal communication habits of two different users 
independently chosen from the same frequency classes, as in¬ 
troduced in Fig.|3] The coefficient covers 3 potential relations 
in the communication interactions, e.g. full and RT spike 
trains, red (left) triangles, full and @, blue (right) triangles, 
and finally RT and @, black and green diamonds. These 3 
coefficients are calculated for who (a, b) and whom users (c, 
d), separately. 6 coefficients in total prove that the temporal 
patterns present high consistency in each average frequency 
classes, the activity (an) and the popularity {pu)- In (b, d), 
the corresponding coefficients are described with the sensi¬ 
tivity of the frequency classes since the average frequency in 
the class of RT is so similar, but not exactly equal to that 
of @. The colored areas are as defined in Fig. [5] and charac¬ 
terize the three main regions of the temporal patterns of the 
individual user spike trains, e.g. bursts (orange), irregular 
random (purple), and regular patterns (yellow). 


A. Nomenclature 

• OSM: Online Social Media, 

• Mention a user name in a tweet message, 

• RE: Reply to a tweet or retweet message, 

• RT: Retweet, share a message of other users in 
her/his own tweet blog, 

• WHO: Twitter users starting an interaction via @ 
or RE or RT with any other users, 

• WHOM: Twitter users addressed by the who users 
such that their message is retweeted or user name 
is mentioned in a message by the who users or 
they get a reply from the who users. 

Any relation between who and whom users such as 
the following-follower is not imposed. 


Figure 7. Linear correlations of Lv of the same users: The 
procedure and representation of the coefficients follow the 
same procedure as introduced in Fig. [6l However, we now 
impose the same users in the same frequency classes. Even 
though (a, c) present the agreement in the temporal patterns 
of full and RT spike trains of the same users, with high corre¬ 
lation coefficients in almost all frequency ranges, (b) indicates 
lower consistency between RT and @ spike trains during en¬ 
tire activity {au) and (d) provides a significant result. While 
less temporal coherence is observed between RT and @ spike 
trains in low popularity {pu), the correlation drastically in¬ 
creases with {pu)- 


V. DISCUSSION 

In this paper, our interest is to quantify online user 
communication in Twitter. To reduce the complexity in 
the communication, the data studied here consider only 
a unique subject which users talk about such as the dis¬ 
covery of the Higgs boson on July 4, 2012 within a re¬ 
stricted time window, e.g. 6 days [^. The main aim 
is to extract salient temporal patterns of communication 
in various types of interaction observed in Twitter such 
as retweet (RT), mention (@), and reply (RE). Adopt¬ 
ing the technique so-called local variation Ly originally 
introduced for neuron spike trains (39l - l^ and rece ntly 
has applied to hashtag spike trains in Twitter [i^, l44l| . 
we perform detail analysis on user communication spike 
trains. Showing strong influences of the frequency of the 
hashtag spike trains on the resultant temporal patterns 
in the earlier work 1 ^ , in parallel we here examine 
the differences in the patterns induced by the frequency 
of the user communication spike trains, fu. 

We investigate user communication spike trains in two 
categorizations, first set of users are the active ones, who 
users, and the other set is composed of the passive users, 
whom users, in the communication. Eor the who users, 
fu simply gives what extend the users contact to the 
whom users and so it is the activity of the who users. 
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au- On the other hand, for the whom users, the gener¬ 
ated spike trains present how often the who users refer 
the messages or the user names of the whom users and 
therefore, fu is the popularity of the whom users, pu- 
Providing comparative statistics on Ly of who and whom 
with increasing ajj and pu, respectively, we observe quite 
distinct temporal behavior of online users. First, we ob¬ 
serve an asymmetry between active and passive inter¬ 
actions, as only the former give rise to hubs, with few 
users attracting a large share of the attention. Moreover, 
who users constantly present bursty patterns, Ly > 1 
for all values of au, whereas whom users demonstrate 
various dynamic behavior patterns, depending on their 
popularity: The least popular users with low pu experi¬ 
ence bursty time series, popular users with moderate and 
high Pu are contacted by temporarily uncorrelated who 
users and so show Poisson random spike trains Ly ss 1, 
and the most popular users with the maximum pu are 
referred regularly in time, e.g. Ly < 1. 

These scenarios are independent of both the position 
of the users, e.g. who or whom, and the preferred inter¬ 
actions, e.g. whether RT or @, suggesting that the fre¬ 
quency of the communication dominates to design social 
dynamic behavior. This conclusion is also supported by 
the high correlation coefficient of Ly on the user pairs 
in the same frequency classes. Furthermore, the linear 
correlation of Ly on the same users reveals interesting 
patterns. There, we observe that only popular users have 
similar dynamic behavior in both RT and which con¬ 
firms that both metrics are complementary to character¬ 
ize the influence of users. 

The analysis could be specified by integrating the com¬ 
munication spike trains with the following-follower rela¬ 
tion in Twitter, and focusing on the who and whom trains 
of connected users. An important concern is the limited 
time period of the data which the collection started 3 
days before the announcement of the discovery and con¬ 
tinued until 3 days after this date. Yet, it has been shown 
that the dynamics of the communication is drasti cally 
different before/after and during the announcement [38l| . 
and this variation could be investigated in our analysis. 
Our study shares the similar aims of the other research 
on online user behavior and the influence of the frequency 


in online platforms such as Flickr, Delicious and Stum- 
bleUpoUjwhich user profiles have been included in the 
analysis [43. This understanding could be also applied to 
our analogy with considering further details in the data. 


A. Data Sharing 

The full data studied in this paper has open access [H, 

Eg. 

DISCLOSURE/CONFLICT-OF-INTEREST 

STATEMENT 

The authors declare that the research was conducted in 
the absence of any commercial or financial relationships 
that could be construed as a potential conflict of interest. 


AUTHOR CONTRIBUTIONS 

Conceived and designed the experiments: CS. Per¬ 
formed the experiments: CS. Analyzed the data: CS. 
Contributed reagents/materials/analysis tools: RL CS. 
Wrote the paper: CS RL. 


ACKNOWLEDGMENTS 

C. Sanh acknowledges supports from the European 
Union 7th Framework OptimizR Project and FNRS (le 
Fonds de la Recherche Scientifique, Wallonie, Belgium). 
This paper presents research results of the Belgian Net¬ 
work DYSCO (Dynamical Systems, Control, and Opti¬ 
mization), funded by the Interuniversity Attraction Poles 
Programme, initiated by the Belgian State, Science Pol¬ 
icy Office. 

Funding. The EU 7th Framework OptimizR Project: 
48909A2 CE OPTIMIZR (Grant holder: RL, Eunding 
receiver: CS - http://optimizr.eu/) and E.N.R.S MIS 
E4527.I2 48888F3 (Grant holder: RL, Funding receiver: 
GS - http://www.fnrs.be/). 


[1] A. Abisheva, V. R. K. Garimella, D. Garcia, and I. We¬ 
ber, Proceedings of the 7th ACM International Confer¬ 
ence on Web Search and Data Mining , 593 (2014). 

[2] S. Wu, J. M. Hofman, W. A. Mason, and D. J. Watts, 
Proceedings of the 20th International Conference on 
World Wide Web , 705 (2011). 

[3] Q. Xuan, M. Gharehyazie, P. Devanbu, and V. Filkov, 
Proceedings of the 2012 International Conference on So¬ 
cial Informatics , 78 (2012). 

[4] M. D. Choudhury, Gamon, Counts, and Horvitz, Pro¬ 
ceedings of the Seventh International AAAI Conference 
on Weblogs and Social Media , 128 (2013). 


[5] R. Lambiotte and M. Kosinski, Proceedings of the IEEE 
102 , 1934 (2014). 

[6] J. Mathiesen, L. Angheluta, P. T. H. Ahlgren, and M. H. 
Jensen, Proceedings of the National Academy of Sciences 
110 , 17259 (2013). 

[7] E. Bakshy, J. M. Hofman, W. A. Mason, and D. J. Watts, 
Proceedings of the Fourth ACM International Conference 
on Web Search and Data Mining , 65 (2011). 

[8] S. Gonzalez-Bailon, J. Borge-Holthoefer, and Y. Moreno, 
American Behavioral Scientist 57, 943 (2013). 

[9] E. Bakshy, B. Karrer, and L. A. Adamic, Proceedings 
of the 10th ACM Conference on Electronic Commerce , 




325 (2009). 

[10] Y.-R. Lin, J. P. Bagrow, and D. Lazer, Proceedings of 
the Fifth International AAAI Conference on Weblogs and 
Social Media , 193 (2011). 

[11] K. H. Kwon, M. A. Stefanone, and G. A. Barnett, Amer¬ 
ican Behavioral Scientist 58 , 1345 (2014). 

[12] J. Zhang, J. Tang, J. Li, Y. Liu, and C. Xing, ACM 
Trans. Knowl. Discov. Data 9, 25:1 (2015). 

[13] K. Lerman, X. Yan, and X.-Z. Wu, ArXiv e-prints 
(2015), arXiv:1506.03022 

[14] A.-L. Barabasi, Nature 435 , 207 (2005). 

[15] K.-I. Goh and A.-L. Barabasi, EPL (Europhysics Letters) 
81 , 48002 (2008). 

[16] G. Miritello, R. Lara, and E. Moro, Temporal Networks 
, 175 (2013). 

[17] M. Formentin, A. Lovison, A. Maritan, and G. Zanzotto, 
Phys. Rev. E 90, 012817 (2014). 

[18] P. Wang, T. Lei, C. H. Yeung, and B.-H. Wang, EPL 
(Europhysics Letters) 94, 18005 (2011). 

[19] T. Zhou, Z.-D. Zhao, Z. Yang, and C. Zhou, EPL (Eu¬ 
rophysics Letters) 97 , 18006 (2012). 

[20] E. R. Colman and D. Vukadinovic Greetham, Phys. Rev. 
E 92 , 012817 (2015). 

[21] M. Karsai, K. Kaski, A.-L. Barabasi, and J. Kertesz, Sci. 
Rep. 2 , 397 (2012). 

[22] H.-H. Jo, J. 1. Perotti, K. Kaski, and J. Kertesz, 
Phys. Rev. E 92 , 022814 (2015) 

[23] H.-H. Jo, M. Karsai, J. Kertesz, and K. Kaski, New 
Journal of Physics 14 , 013055 (2012). 

[24] J. Kim, D. Lee, and B. Kahng, PLoS ONE 8, e58292 
(2013). 

[25] J. L. Iribarren and E. Moro, Phys. Rev. Lett. 103, 038702 
(2009). 

[26] R. Lambiotte, L. Tabourier, and J.-C. Delvenne, The 
European Physical Journal B 86 (2013). 

[27] C. Wang and B. A. Huberman, Sci. Rep. 2, 633 (2012). 

[28] M. Formentin, A. Lovison, A. Maritan, and G. Zanzotto, 
ArXiv e-prints (2014), arXiv: 1405.5726 [physics.soc-ph] 

[29] E. Ferrara, R. Interdonato, and A. Tagarelli, Proceedings 
of the 25th ACM Conference on Hypertext and Social 
Media , 24 (2014). 

[30] J. Ratkiewicz, S. Fortunato, A. Flammini, F. Menczer, 
and A. Vespignani, Phys. Rev. Lett. 105 , 158701 (2010). 

[31] M. Coscia, International AAAI Conference on Weblogs 
and Social Media (2013). 

[32] S. A. Myers and J. Leskovec, Proceedings of the 23rd In¬ 
ternational Conference on World Wide Web , 913 (2014). 

[33] Algorithms and Models for the Web Graph 8882, 132 
(2014). 

[34] Q. Zhao, M. A. Erdogdu, H. Y. He, A. Rajaraman, and 


J. Leskovec, KDD15 (2015). 

[35] D. Boyd, S. Colder, and G. Lotan, System Sciences 
(HICSS), 2010 43rd Hawaii International Conference on 
, 1 ( 2010 ). 

[36] T. R. Zaman, R. Herbrich, J. V. Gael, and D. Stern, 
Computational Social Science and the Wisdom of Crowds 
Workshop (colocated with NIPS 2010) (2010). 

[37] D. V. Greetham and Ward, Proceedings of the 2nd In¬ 
ternational Workshop on Dynamic Networks and Knowl¬ 
edge Discovery , 73 (2014). 

[38] M. D. Domenico, A. Lima, P. Mougel, and M. Musolesi, 
Sci. Rep. 3, 2980 (2013). 

[39] S. Shinomoto, K. Shima, and J. Tanji, Neural Comput. 
15, 2823 (2003). 

[40] K. Miura, M. Okada, and S. ichi Amari, Neural Comput. 
18 , 2359 (2006). 

[41] T. Omi and S. Shinomoto, Neural Comput. 23, 3125 

( 2011 ). 

[42] S. Shinomoto, H. Kim, T. Shimokawa, N. Matsuno, 
S. Funahashi, K. Shima, 1. Fujita, H. Tamura, T. Doi, 

K. Kawano, N. Inaba, K. Fukushima, S. Kurkin, K. Ku- 
rata, M. Taira, K.-I. Tsutsui, H. Komatsu, T. Ogawa, 
K. Koida, J. Tanji, and K. Toyama, PLoS Comput Biol 
, el000433 (2009). 

[43] C. Sanli and R. Lambiotte, PLoS ONE 10, e0131704 
(2015). 

[44] C. Sanli and R. Lambiotte, Modeling and Mining Tem¬ 
poral Interactions: Papers from the 2015 ICWSM Work¬ 
shop , 8 (2015). 

[45] J. Leskovec and A. Krevl, SNAP Datasets: Stanford 
Large Network Dataset Collection (2014). 

[46] J. Saramaki, E. A. Leicht, E. Lopez, S. G. B. Roberts, 
F. Reed-Tsochas, and R. 1. M. Dunbar, Proceedings of 
the National Academy of Sciences 111 , 942 (2014). 

[47] P. d. Meo, E. Ferrara, F. Abel, L. Aroyo, and G.-J. 
Houben, ACM Trans. Intell. Syst. Technol. 5, 14 (2014). 

[48] B. A. Huberman, D. M. Romero, and F. Wu, First Mon¬ 
day 14 (2009). 

[49] B. Goncalves, N. Perra, and A. Vespignani, 
PLoS ONE 6, e22656 (2011) 

[50] L. Weng, A. Flammini, A. Vespignani, and F. Menczer, 
Sci. Rep. 2 , 335 (2012). 

[51] J. P. Gleeson, J. A. Ward, K. P. O’Sullivan, and W. T. 
Lee, Phys. Rev. Lett. 112, 048701 (2014). 

[52] U. Cetin and H. O. Bingol, Phys. Rev. E 90, 032801 
(2014). 

[53] J. Borge-Holthoefer, R. A. Banos, S. Gonzalez-Bailon, 
and Y. Moreno, Journal of Complex Networks (2013). 

[54] H.-H. Jo, R. K. Pan, and K. Kaski, Phys. Rev. E 85 , 
066101 (2012). 


